Self-Doping Induced Orbital-Selective Mott Transition in Hg 2 Ru 2 07 
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Pyrochlore oxides are fascinating systems where strong, multi-orbital correlations in concert with 
geometrical frustration give rise to unanticipated physical properties. The detailed mechanism of 
the insulator-metal transitions (IMT) underpinning these phenomena is, however, ill-understood in 
general. Motivated thereby, we study the IMT in the pyrochlore Hg2Ru2C>7 using LDA+DMFT. 
In contrast to the well-known examples of Mott transitions in TMO, we show that, in the negative 
charge-transfer situation characteristic of Hg2Ru207, self-doping plays a crucial role in the emer- 
gence of an orbital-selective IMT. We argue that this mechanism has broader relevance to other 
correlated pyrochlore oxides. 

PACS numbers: 71.28+d,71.30+h,72.10-d 
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I. INTRODUCTION 

The Mott-Hubbard insulator-metal transition (IMT) 
is by now recognized to play a central role in our un- 
derstanding of d- and /-band compounds. 1 Understand- 
ing the complex interplay between strong, multi-orbital 
(MO) electronic correlations, structural distortions, and 
strongly anisotropic, orbital-dependent hopping holds 
the key to a consistent theoretical description of their 
unique responses. Adding geometric frustration to the 
above results in a truly formidable problem. In geo- 
metrically frustrated systems, the exponentially large de- 
generacy of classical ordered states inhibits emergence 
of conventional order, permitting new, complex ordered 
ground states to arise^ In real systems, structural ef- 
fects may partially remove this huge degeneracy, making 
the problem (counter-intuitively) somewhat simpler to 
solved however, in near-undistorted cases, near perfect 
orbital degeneracy, and the consequent strong quantum 
orbital and spin fluctuations in a highly degenerate sys- 
tem underpin their physical behavior. 

The recently discovered pyrochlore system, 
Hg2Ru2C>7, 4 is a particularly interesting case in 
this context. As temperature (T) is reduced, the 
anomalous (see below) non-Fermi liquid (nFL) metallic 
state becomes unstable, via a first-order Mott transition, 
to an antiferromagnetic Mott-Hubbard insulator (AFI). 5 
External pressure (p) melts the AFI beyond a critical 
p c = 6.0 GPa, resulting in a low-T correlated FL 
behavior for T < T* = 13K. At ambient pressure, the 
IMT transition is accompanied by lowering of lattice 
symmetry from (high-T) cubic to (low-T) a lower, 
hitherto precisely unknown, type: depending upon its 
precise type, either AF or dimer order may be possible. 6 

The high-T > T* state in Hg2Ru207 is an anomalous 
nFL: the dc resistivity, p(T) ~ T for T > T MI = 108 K± 
at ambient pressure, and deviates from the FL form for 
T > T* beyond p c . The uniform spin susceptibility is 
Curie- Weiss like for T > Tmi, indicating a strong local 
moment scattering regime. The nFL character is also 
borne out from the recent photoemission (PES) data, 7 



showing anomalously broad PES lineshapes, with no hint 
of FL quasiparticles, in the metallic phase. Given the 
cubic pyrochlore structure for T > Tmi, strong orbital 
(from t2 g orbital degeneracy) and spin fluctuations are 
implied. How might strong scattering processes involv- 
ing these fluctuations produce the observed nFL metal? 
What drives the AFI as T is lowered? A correlated elec- 
tronic structure study which can illuminate these issues 
does not, to our best knowledge, exist. In this work we 
study precisely these issues in part in Hg2Ru2C>7 using 
the LDA+DMFT method.— We focus on the mechanism 
of the T-driven IM transition, and leave the issue of the 
low-dimensional AF with a spin gap for future consider- 
ation. 



II. MODEL AND SOLUTION 

Starting with the high-T cubic Fd3m structure found 
by Klein et alA local density approximation (LDA) band 
structure calculations were performed using a (scalar 
and fully-relativistic) full-potential local-orbital scheme 
(FPLO) 9 and a linear muffin-tin orbitals (LMTO) scheme 
in the atomic sphere approximation. 10 In Fig. Q] we 
display our FPLO-LDA 11 results for the one-particle 
density-of-states (DOS). Clearly, the major contribution 
to the DOS comes from Ru(Ad), 0(l)(2p) orbitals, but 
the 0(2) (2p) orbitals also have noticeable spectral weight 
at the Fermi level (Ep). As seen in Fig.[IJ the influence 
of the spin-orbit coupling (SOC) in the cubic phase to 
the tig states near Ep is negligible. Further, we ob- 
serve strong hybridization between Ru(4d) — 0(1) (2p) 
orbitals, and between 0(2)(2p) — Hg(6s) orbitals, but 
weak mixing between these two sets. These important 
findings will be exploited below to study the physics 
of Hg2Ru207Using a multi-orbital Hubbard model in- 
volving only the d-band sector. At one-particle level, 
the corresponding model Hamiltonian reads, Hband = 
J2k,a,a tkadi aa d kaiJ + £»,o,6=i,2 A 6 (n? 6 - nf a ), where a 
labels the diagonalized combination of (e' gl , e' g2 , a-ig) 12 
the t 2g orbitals, (b = 1,2) are the pd charge trans- 
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for (CT) terms involving two Ru — 0(1), (b = 1) and 
Ru — 0(2), (b — 2) channels. Clearly, neither an AFI 
Mott-Hubbard insulator nor a nFL metal can be expected 
at this level, this requiring a reliable treatment of strong, 
c?-shell electronic correlations. This part reads, 
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with a, a' — e' gl ,e' g2 ,ai g . Given the larger spatial ex- 
tent of 4d-orbital vis-a-vis their 3d counterpart, we also 
include a Madelung term, H M = U pd Y, <it j >iaib nf ta n^ b 
in our calculations (see below). In Tl 2 Mn 2 0y, the IMT 
(from a paramagnetic insulator (PI) to a ferromagnetic 
metal (FMM)) is seemingly driven by the Tl(6s) states 
crossing Ep across T cf H for instance. This can occur via 
a T-dependent CT from the TM d-states to the 0(2)-Tl 
hybridized states. In view of the generic relevance of self- 
doping in TM 4<i pyrochlores 1 13 i 14 this term is an essential 
part. In the low-T phase, this will involve CT processes 
involving the 0(2p) — Hg(6s) channel (second term of 
Hband) across the IMT, as we describe below. 
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FIG. 1: (Color online) LDA band structure for cubic 
Hg2Ru207. 0(1) (0(2)) denotes Oxygen ions nearest to Ru 
(Hg). 

The daa above should be understood as appropriate, 
Ru0 6 cell-centered, combinations of the Ru-4d and O- 
2p orbitals, computed within LDA. Similarly, the b = 2 
channel is to be understood as a band of hybridized O- 
2p and Hg-6s orbital states. Because of the negligible 
one-particle mixing (hybridization) between the b = 1,2 
channels, we approximate the full problem of three Ru-c£ 
bands coupled to the b — 2 band channel by replacing the 
latter by a "reservoir" , whose only function is to simulate 
the self-doping process arising from the negative charge 
transfer situation that obtains in Hg2Ru207. Of course, 
this is an approximation. It is, however, a good one: the 



b = 2 band channel has appreciably smaller DOS around 
Ep in the LDA results. In a DMFT-like approxima- 
tion, with negligible one-particle hybridization between 
the 6 = 1,2 bands, the inter site Madelung term will push 
this small spectral weight away from Ep, to lower- and 
higher energy (notice that the 0-2p bands will be split 
by U p dz{nd), where z is the co-ordination number of the 
lattice, and hence quite large). This will already occur 
at the level of LDA+Hartree approximation. We then 
expect that the correlated spectral function will be dom- 
inated by the c£-bands over an appreciable range about 
E F . Correlation effects in Hg 2 Ru 2 7 via LDA+DMFT 
are studied below, subject to this caveat; we will show 
that this is indeed a good approximation a posteriori, in 
the sense that our LDA+DMFT results with the above 
caveat show very good quantitative agreement with one- 
particle spectroscopy and key thermodynamic and trans- 
port data in Hg 2 Ru 2 07. 

The full many-body Hamiltonian reads H = Hb an d + 
Hint + Hm- We solve this model within MO dynamical 
mean field theory (MO-DMFT) developed and used for 
a range of TMO with good success^ We use the MO- 
iterated perturbation theory (IPT) as an impurity solver 
in the DMFT selfconsistent procedure. Though not 
numerically exact (like QMC, NRG, D-DMRG), it has 
many advantages: it is valid for T — 0, where QMC can- 
not be used. NRG/DDMRG are extremely prohibitive 
for three orbital models, even without electronic struc- 
tural inputs at LDA level. As shown in earlier work, 15 
DMFT(MO-IPT) generically gives very good semiquanti- 
tative agreement with PES/X-ray absorption (XAS) ex- 
periments for TMOsJ^ For Hg 2 Ru 2 7 , DMFT(MO- 
IPT) has to be extended to treat the second CT chan- 
nel described above. Given the complexity of the 
problem, we choose the following strategy to solve H 
above: For the high-T phase, we (i) solve Hband + H int 
within DMFT (MO-IPT) version used earlier, for tech- 
nical details see Ref. |l5|. (ii) As indicated by ex- 
periment on the related Tl 2 Ru 2 0y system^ we study 
the effect of Ru(4d) - 0(2)(2p) - Hg(6s) CT processes 
by incorporating this CT channel selfconsistently into 
the LDA-t-DMFT(MO-IPT) procedure: given the small 
Ru(4d) - 0(2)(2p) hybridization, the 0{2){2p) - Hg(6s) 
channel acts like a scattering (non-hybridizing) channel 
for the Ru — i 2g bands in the impurity model of MO- 
DMFT. This enables us to treat this extra channel by 
combining the MO-IPT solution for (i) with the exact 
DMFT solution of a Falicov-Kimball model (FKM) for 
(ii)r± in a selfconsistent way. 



III. LDA+DMFT RESULTS AND DISCUSSION 

We now describe our LDA+DMFT results. We start 
with the (cubic) pyrochlore structure at high-T, with the 
corresponding LDA DOS as the input for the MO-DMFT 
calculation with total Ad occupation, n t — 3. Further, 
we work in the (LMTO)^ representation in which the 
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one-particle density matrix is diagonal in the t 2g orbital 
index, so that G^(k,u) = 5 a pGal{k,uj). We choose 
U = 5.5 eV, J H = 1.0 eV, and U' ~ (U - 2J„) = 3.5 eV 
for the Ru — Ad shell, along with U p d — 1.5 eV, as repre- 
sentative values for Hg2Ru2 07. We believe that U, Jh 
do not vary much for 4d-TMO pyrochlores, and, in fact, 
U = 5.0 eV was found for the CT insulator Cs 2 AgF 4 ^ 




a(eV) 



FIG. 2: (Color online) t 29 -resolved LDA+DMFT densities of 
states for Hg2Ru207in the high-T nFL phase, for U — 5.5 eV 
and Jh = 1 eV. Notice the orbital selective nature of the 
metal, as well as the nFL character in the orbital resolved 
self-energies (see inset). 



A. Metallic Phase 

In Fig. [21 we show the correlated many-particle spec- 
tral function for the metallic phase of Hg2Ru207. The 
dynamical spectral weight transfer (SWT) over large en- 
ergy scales, characteristic of strong, local correlations is 
explicitly manifest. More interestingly, the metal has an 
orbital selective (OS) character: the a\ g DOS is almost 
"Mott" localized, even as the e' g DOS develops a pre- 
cursor of a low-energy pseudogap, characteristic of an 
incoherent metal behavior. This is further corroborated 
by examining the orbital dependent self-energies (T,(u>)), 
see inset in Fig. [5] Clearly, the e' g (imaginary part of) 
E(w) shows quasi-linear frequency dependence near the 
Fermi energy (\u> — Ep | <0.3 eV), while ImS ai (u>) shows 
quadratic behavior. This describes a nFL metal, with a 
linear-in-T quasiparticle damping rate. Within DMFT, 
this is also the transport relaxation rate^ since vertex 
corrections drop out in the computation of the conduc- 
tivities in this limit. The dc resistivity is then given by 
p dc {T) ~ (m* /ne 2 )lm£ e , g (w = T) ~ AT, in accord with 
the linear-in-T resistivity observed experimentally in the 
"high-T" metallic phase. Moreover, the selective localiza- 



tion seen in the a% g orbital DOS gives rise to unquenchcd 
local momemts co-existing with "itinerant" (but incoher- 
ent as derived above) e' g carriers naturally gives a Curie- 
Weiss form of the spin susceptibility, which is also ob- 
served right up to the IMT. We have also estimated the 
7-co-efficient of the metallic specific heat from the real 
part of T, a (uj) (not shown) as j/^lda = 4.25: this seems 
to be in the range estimated in Ref 0. Actually, the 
noticeable T dependence of 7 above TmA is additional 
evidence of disordered local moments in the "bad" metal, 
and further supports our picture. 
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FIG. 3: (Color online) Theoretical PES spectra 
(LDA+DMFT DOS convoluted with instrumental res- 
olution) in the nFL metallic and insulating phases of 
Hg 2 Ru 2 7 . In the LDA+DMFT calculation, the total t 2g 
electron density changes from < n t >= 3.0 (metallic) to 
< n t >— 2.6 (insulating). The theoretical result shows very 
good agreement with experimental PES result in the metallic 
phase, taken from Ref. |7| up to 2 eV binding energy. The 
inset shows the total spectral function. 

Using the LDA+DMFT result, we also compare 
(Fig. [3]) our computed PES lineshape, with very recent 
work from the RIKEN-SPring 8 group^ Given our re- 
striction to the t 2g - Ru(4d) bands in the MO-DMFT 
(Hg(d),0(2p) bands will start contributing at higher 
binding energies, as seen from LDA), good quantitative 
agreement with experiment is evident up to —2.0 eV, 
lending strong support to our theoretical work. Addi- 
tionally, we predict that an intense peak will be seen 
around 0.4 — 0.5 eV in X-ray absorption (XAS) studies, 
as in Tl 2 Ru 2 7 ^ 

In fact, the LDA+DMFT spectral functions show that, 
at low energy the full MO problem is mapped onto an ef- 
fective FKM-like models, since the bad-metal is of the 
orbital-selective type. There, the problem in the local 
limit corresponds to itinerant (but incoherent) e' g car- 
riers scattering off (Mott) localized a± g electronic states. 
The resulting problem is precisely the "X-ray edge" prob- 
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lem in DMFTj^ 2 - We now understand the structure of 
the self-energies, and the nFL behavior, as an interest- 
ing manifestation of the Anderson OC caused by this 
X-ray-edge mapping. Building upon this understanding, 
using DMFT, we predict that, for T > T MI : (a) the op- 
tical conductivity will show a low-energy pseudogapped 
form, characteristic of an incoherent (nFL) metal. Po- 
larised optical studies should indicate the OS-character 
of the nFL state, and, (b) given the local version of the 
Shastry-Shraiman relation^ the electronic Raman line- 
shape should show a continuum response, cut off by a 
pseudogap feature at low energy. 



B. Insulating Phase 

Now, we turn to a description of the insulating phase. 
A naive search for the instability of the nFL to a Mott- 
Hubbard insulating state, where U, U' were increased to 
unphysical values, nevertheless turned out to be unsuc- 
cessful. This is a clear indication of involvement of ad- 
ditional, electronic-cum-structural effects in driving the 
IMT. The observed structural change across Tmi^ sup- 
ports this reasoning. From the LDA orbital assignment, 
it is clear that an additional structural change necessar- 
ily involves partial occupation of the two-fold degener- 
ate e' g orbitals. Starting from the nFL-metal derived 
above, this can only occur via a T-dependent change 
in the "self-doping" process.— As in TI2RU2O7, this In- 
dependence could be provided by electronic coupling to 
the Ru — O stretching phonon mode, which is experimen- 
tally observed to split below Tmi-^ This would imply a 
structural change setting in below Tmi, whose precise 
nature is hitherto unknown. 

Motivated by this observation, we argue that this re- 
sulting change in the Ru(4d) — > 0(2) (2p) CT leads to a 
partial Ad occupation, induces orbital order, and lifts the 
e' g degeneracy via a structural distortion. Within DMFT, 
this will reduce the orbital-dependent hoppings, driving 
large SWT from low- to high energy, and stabilize the 
second, Mott-Hubbard insulating, solution of the DMFT 
equations. This is indeed seen in the LDA+DMFT cal- 
culation, as we show below. 

Based on this reasoning, we explore the Mott-Hubbard 
insulating phase of Hg2Ru207by searching numerically 
for the instability of the first (metallic) solution found 
above to the second, insulating solution of the DMFT 
equations in the quantum paramagnetic phase. The 
DMFT(MO-IPT) equations are solved with U,U',U pd , 
and nt = 2.6. We vary n t in trial steps, and look for a 

(c) 

critical n t which stabilizes the PI solution of the DMFT 
equations. As discussed above, and indeed seen in the 
inset of Fig. partial occupation of the Ad — e' g or- 
bitals leads to removal of the orbital degeneracy, and the 
corresponding modification of orbital orientation gives 
rise to a reduced intersite one-electron overlap. Within 
MO-DMFT, this triggers the Mott-Hubbard insulating 
state via large SWT on a scale of 5.0 eV, as seen in 



Fig. H Hence, the IMT is an OS Mott transition. Im- 
portantly, however, notice that the self-doping process 
leading to fractional d-orbital occupation is a crucial in- 
gredient. Thus, in contrast to other OS cases^ the IMT 
in Hg2Ru2 0y is driven by the self-doping in the negative 
CT situation, A 2 = {e p2 - e d )=(-2.8eV+1.3eV)=-1.5eV. 
Along with the microscopic elucidation of the nFL be- 
havior, the agreement with PES in the nFL metal phase 
up to lu ~ —2.0 eV (Fig. [3]) constitutes strong evidence 
in favor of our mechanism for the IMT in Hg^R^Oy. 
Quantitative comparison with PES in the low-T phase 
requires an extension of our approach to include short- 
ranged (intersite) orbital and spin correlations charac- 
teristic of pyrochlores. This requires a cluster-DMFT 
analysis, presently a forbidding prospect. Nevertheless, 
observation of a Curie- Weiss spin susceptibility right up 
to the IMT,— i2a along with the excellent agreement with 
PES in the nFL, justifies the use of DMFT to describe 
the IMT. 

We emphasize that this is a new picture for the IMT in 
correlated systems. In contrast to the early TMO, which 
are Mott-Hubbard insulators {Udd > A in the Zaanen- 
Sawatzky-Allen scheme)^ the Ru(4d) - 0(1) (2p) CT 
channel, along with the Madelung term, are important 
in Hg2Ru207. In the DMFT context, the importance 
of the CT energy in the late-TMO is recognized. 26 Our 
work is the first of its kind showing its relevance for the 
IMT in a MO pyrochlore system. Further, it is likely 
to be more broadly generic to pyrochlore TMOs; recall 
that, in Tl 2 Mn 2 07, the PI-FMMi^ transition is driven 
by the shift of the Tl(6s) band through Ep across T c . 
This is readily rationalized in our picture, in terms of 
the "switching on" the Mn(3d) - 0(2) (2p) - TZ(6s) CT 
channel across T c in Tl 2 Mn 2 07F^ 

The correlation between the CT process detailed above 
and the IMT is also visible in a whole family of id py- 
rochlores, A 2 Ru 2 7 , with A = Pb, Bi, Y£- The low- 
T magnetic ordering in the Mott-Hubbard insulating 
phase(s) may, however, be quite different, being sensi- 
tively dependent on the nature of the structural change 
across the IMT, as well as on spin state. For example, 
the spin 5 = 1 system, TI2RU2O7, shows a very similar 
Mott transition; however, the low-T phase is found to be 
a Haldane spin chahii^Z In Hg2Ru207, the half-integer 
spin 5 = 3/2 rules out the Haldane analogy. If the low- 
T magnetic structure corresponds to having Ru chains, 
as in TI2RU2O7, one would have an AF ground state 
with gapless spin excitations. Observation of the spin- 
gap in the uniform spin susceptibility in Hg2Ru207 may 
therefore point to the relevance of the spin-orbit cou- 
pling in the insulating phase: this will induce Ising-like 
anisotropy in an 5 = 3/2 Heisenberg chain and generate 
a gap to spin excitations^ More experimental work is 
also called for to pinpoint the specific factors affecting 
this issue. Given the structural distortion necessarily ac- 
companying the IMT caused by lifting of the e' g degener- 
acy in Hg 2 Ru 2 07, more detailed theoretical work awaits 
more precise characterization of the low-T structure of 



5 



Hg 2 Ru207- We plan to address the issue of magnetic or- 
der and its associated excitation spectrum in the low-T 
(insulating) phase of Hg 2 Ru207 in a future work. 

IV. CONCLUSION 

In conclusion, we have performed LDA+DMFT calcu- 
lations to demonstrate the role of multi-orbital electron- 
electron interactions on the electronic structure of a re- 
cently discovered pyrochlore system, Hg2Ru207. This 
system exhibits a first-order temperature dependent 
IMT. The high-T metallic phase is shown to be an orbital 
selective non-Fermi liquid, signalled by a quasi-linear fre- 
quency dependence of the imaginary part of the corre- 
lated self-energy and the absence of quasi-particle peaks 
in the orbital-selective spectral functions at very low 
frequencies. The LDA+DMFT spectral function shows 
good quantitative agreement with recent photoemission 
(PES) data, showing anomalously broad PES lineshapes, 



as well as with the linear-in-T resistivity and mass en- 
hancement in the metallic phase. In stark contrast to 
the better known examples 1 * of Mott transitions in 3d- 
transition metal compounds, our results imply a new 
mechanism for MIT in the 4d pyrochlore-based electron 
systems (like Hg 2 Ru 2 07). Namely, a negative charge- 
transfer associated self-doping drives an orbital-selective 
MIT at low temperatures via the Mott-Hubbard route. 
Our study should be more generally applicable to MO 
pyrochlore systems, and, in particular, to Tl2Mn 2 07, ex- 
hibiting Mott-Hubbard transition o 14 ' 21 ' 24 ' 29 as functions 
of suitable "tuning parameters" . 
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